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Renormalization Group Technique for 
Quasi-one-dimensional Interacting Fermion 
Systems at Finite Temperature 

C. Bourbonnais, B. Guay and R. Wortis 

ABSTRACT We review some aspects of the renormalization group method for interacting fermions. Special emphasis 
is placed on the application of scaling theory to quasi-one-dimensional systems at non zero temperature. We begin by 
introducing the scaling ansatz for purely one-dimensional fermion systems and its extension when interchain coupling 
and dimensionality crossovers are present at finite temperature. Next, we review the application of the renormalization 
group technique to the one- dimensional electron gas model and clarify some peculiarities of the method at the two- 
loop level. The influence of interchain coupling is then included and results for the crossover phenomenology and 
the multiplicity of characteristic energy scales are summarized. The emergence of the Kohn-Luttinger mechanism in 
quasi-one-dimensional electronic structures is discussed for both superconducting and density-wave channels. 
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1 Introduction 



Scaling ideas have exerted a far-reaching influence on our understanding of complex many-body systems. Their 
use in the study of critical phenomena has offered a broad and fertile field of work at the heart of which is the 
Wilson formulation of the renormalization group method. |jl|, ^ In the Wilson view, fluctuations of the order 
parameter at all length scales up to the correlation length are the key ingredient in explaining the existence of 
power law singularities that govern critical properties in accordance with scaling laws. Equally powerful is the 
extension of scaling concepts to anisotropic systems, namely when small anisotropic parameters are magnified as 
a result of their coupling to singular fluctuations. Intermediate length scales thus emerge giving rise to changes or 
crossovers in the critical behavior.]^, ^, |5) The horizon of applications of this methodology was further widened 
when scaling was applied to quantum critical systems for which quantum mechanics governs fluctuations of the 
order parameter at the critical point. H, 0] 

There has been a parallel expansion of the use of scaling tools in the description of many-fermion systems. 
This was especially true for the Kondo impurity and the one-dimensional fermion gas problems. In the latter, 
the reduction of spatial dimension is well known to enlarge the range of fluctuations which have a peculiar 
impact on the properties of the system. [||, ^ In this context the scaling theory commonly termed multiplicative 
renormalization group, 1 10, ll[ contributed significantly to the completion of a coherent microscopic picture of the 
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ID fermion gas system as a paradigm of a non Fermi — Luttinger — liquid state. However, this scaling theory 
rests for the most part upon the logarithmic structure of infrared singularities that compose ID perturbation 
theory and as such departs somewhat from the standard Wilson picture. It remains closer to earlier formulations 
of the renormalization group in quantum field theory. Yet, fluctuations that characterize the ID fermion gas do 
show a multiplicity of length scales stretching from the shortest distance, of the order of the lattice constant, 
up to the quantum coherence length or the de Broglie wavelength of fermions. Early attempts to reexamine 
the properties of low dimensional fermion gas system along these lines were motivated by the coupled chains 



problem, |12|, |13|, |14|, 15, 16| which finds concrete applications in the physics of quasi-one-dimensional organic 



conductors. p/A, p_8| The description of interacting fermions when strong spatial anisotropy is present shares many 
traits with its counterpart in critical phenomena. Q Hence the Wilson method provides a powerful framework 
to study how deviations from perfect ID scaling are introduced by spatial anisotropy. The corresponding length 
scales enter as key components of the notoriously complex description of dimensionality crossovers in which the 
ID state evolves either towards the emergence of a Fermi liquid or the formation of long-range order at finite 
temperature. 

The Wilson renormalization group approach has attracted increasing interest in several areas of the many- 
fermion problem, notably in the description of two-dimensional systems in connection with high-Tc materials, |20| 
p2| Fermi liquids, [p3| p^ , psf and ladder systems, p^, p7| , p8| Within the bounds of this review, it would there- 
fore be impossible to give anything more than a cursory account of the expanding literature in this field |2^, ^ |3l[ . 
With the goal of keeping this paper self-contained for the non-specialist, however, our interest will be more selec- 
tive and will focus on the application of the renormalization group to the quasi-one-dimensional fermion gas at 
finite temperature. [3^ Although many aspects of this problem have already been discussed in detail in previous 
reviews, |l^ it is useful to revisit some issues while in addition discussing features of the Wilson method that 
deserve closer examination but which have as yet never received a detailed investigation. This is especially true in 
regard to the formulation of the method itself, in particular the way successive integrations over fermion states in 
outer momentum shells (the Kadanoff transformation) is carried out when high order calculations are performed. 
We close this review by a detailed discussion of the dual nature of the Kohn-Luttinger mechanism|33 from which 
interchain pairing correlations find their origin in quasi-one-dimensional systems.]^, |37|] With the aid 

of the renormalization group, |3^, ^ we show how this mechanism may yield instabilities of the metallic state 
towards either unconventional superconducting or density-wave order. 

In section II, we introduce scaling notions for fermions in low dimension from a phenomenological standpoint 
and predictions of the scaling hypothesis are compared to the solution in the non-interacting case in section III. 
The formulation of the renormalization group in the classical Wilson scheme is given in section IV and explicit 
calculations are carried out up to the two-loop level. In section V, the influence of interchain coupling is discussed 
and we go through various possibilities of dimensionality crossovers in the quasi-one-dimensional Fermi gas. In 
section VI, we study how the Kohn-Luttinger mechanism for interchain pairing correlation emerges from the 
renormalization group flow and under what conditions it can lead to long-range order. We conclude in section 
VII. 



2 Scaling ansatz for fermions 

Before embarking on the program sketched above, it is first useful to discuss on purely phenomenological grounds 
how scaling notions apply to ID and quasi-lD fermion systems. The predictions of the scaling ansatz are explicitly 
checked in the trivial but instructive case of the non-interacting Fermi gas. 



2.1 One dimension 

Let us look at how scaling ideas apply to the description of purely one-dimensional fermionic systems at low 
temperature. The temperature T is the parameter that controls the quantum delocalization or the spatial co- 
herence of each particle in one dimension. This appears as a characteristic length scale, ^, which is precisely 
the de Broglie wavelength. The temperature dependence of ^ is readily determined by considering the effect of 
thermal excitations on fermion states within an energy shell Se{k) ~ T [ks = l,'h= 1), where e{k) is the energy 
spectrum of fermions evaluated with respect to the Fermi level. In the vicinity of the Fermi points ±kp, one has 
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Se{k) « vfH k I — /cf), where vp is the Fermi velocity, and the coherence length can be written 

i^vp/T^T-'', (1.1) 

where the exponent v is then equal to unity. As T — > 0, ^ goes to infinity and the system develops long-range 
quantum coherence. The 'time' r required for quantum delocalization up to ^ is therefore ~ 1/T which simply 
implies that 

r^r. (1.2) 

where z = 1 is the dynamical exponent of the fermion field. Since 'time' and distance must be considered on the 
same footing, both v and z can be taken as independent of interaction. 
Let us now consider the Matsubara-Fourier transform 

G{k,LUn)^ / e-*'=" + ^"""G(a;,r)dxdr, (1.3) 
Jo Jo 

where 

G(x,t) =-(T,^,(x,T)^t (0,0)) (1.4) 

is the one-particle time-ordered correlation function expressed as a statistical average over fermion fields ipa^ of 
spin a. Here LUn = {2n+ 1)t:T corresponds to the Matsubara frequencies. For non interacting fermions and T = 0, 
G{k) '--^ (I fc I — kp)^^ has a simple pole singularity at |fc| kp and = 0, indicating that single particle 
excitations are well defined. However, when fermions interact in one dimension, such a picture is expected to be 
altered. In the scaling picture, this can be illustrated by adding an anomalous power dependence on the wave 
vector near kp^ namely 

G{k) ^ (I fc I - kpr^+r (1.5) 

The absence of a quasi-particle pole introduces 6, the anomalous dimension for the Green's function {9 >{) where 
the equality occurs in the non-interacting limit). Correspondingly, the absence of a Fermi liquid component at 
zero temperature will affect the equal-time decay of coherence at large distance, which takes the form 

G(x) « (1.6) 

where the effective dimension d = 2 has a space and time component. A similar decay in time takes place for 
large r at a: = 0. Looking at the temperature dependence of the Green's function, the scaling hypothesis allows 
us to write 

Gik,iu„)^CT-^g{Sk^,Lu^C), (1-7) 

where Sk = |A;| — kp,j is the thermal exponent of the single particle Green's function, and g{x,y) is a scaling 
function. Consistency between Eqns. (1.7) and ( |l.5[ ) leads to the following relation between the exponents 

7 = (l-0)i/. (1.8) 

Scaling can also be applied to the free energy per unit length, the temperature dependence of which can be 
expressed as 

/ = --InZ w AT^^". (1.9) 
L 

The exponent a is connected to the temperature variation of specific heat 

« A'T-", (1.10) 

where 

a = l-Ji/. (1.11) 
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As long as the ansatz v = I holds for the coherence length, one then always has a = — 1 in one spatial dimension. 
The specific heat is therefore linear in temperature for fermions even though, according to Eqn. ( |l.5| ) , the system 
is no longer a Fermi liquid at low energy. It should be noted that the familiar Sommerfeld result Cy ^ T for a 
Fermi liquid (and Fermi gas) in higher dimensions indicates that the importance put on fermion states near the 
Fermi level one-dimensionalizes the sum over states for the evaluation of internal energy. 

Scaling concepts can also be applied to the critical behavior of correlations that involve pairs of fermions in one 
dimension. Superconducting and 2kp density-wave responses are among the pair susceptibilities that may present 
singularities at zero temperature. Consider the Matsubara-Fourier transform of the susceptibility 

Xiq,u;^) = ft e-''^^+'--^x{x,T)dxdT, (1.12) 



Jo Jo 

where cOm — 2rmTT. The time-ordered pair correlation function is 

x(x,T) = -(T,O(x,r)Ot(0,0)), (1.13) 

in which ^ '4'a'4'p {O^ ^ '^a^is) ^^re the composite operators of the superconducting (density-wave) channel, 
commonly called the Cooper (Peierls) channel. At T = 0, these can show long-range coherence leading to an 
algebraic decay of pair correlations with distance. Following the example of the correlation of the order parameter 
at a critical point, one can write a relation of the form 

X(.)«e---^, (1.14) 

where 77 > is an anomalous exponent that characterizes the decay of pair correlation. Here go = (2kp) is the 
characteristic wave vector of correlations in the Cooper (Peierls) channel. A similar power law variation with time 
T can also be written at small distance. Correspondingly, in Fourier space, one has in the static limit 

xiQ) « ^ , (1-15) 

where q refers to deviations with respect to an ordering wave vector qq in the Cooper or Peierls channel. For 
q = and at finite frequency, a similar expression holds with q being replaced by the frequency. Now at finite 
temperature, the scaling hypothesis allows us to write 

Xix, T) « De"i''''x-^+'^-''X{x/S,), (1.16) 

where X is scaling function (A'(O) = 1 ). Here the correlation length ^ for the pair field is assumed to have the 



same power law dependence on temperature as the coherence length Eqn. (1.1) for both fermions in a pair. In 



Fourier space, one can write for the divergence of the susceptibility in a given channel 

x{q.^.n.T) « DT-^Xiqi^Umn- (1-17) 



Consistency between Eqns. (1.17) and (1.15) in the static case leads to the Fisher relation 

l^{2~ii)y (1.18) 
between the exponents of the static susceptibility, the correlation function and the correlation length. 



2.2 Anisotropic scaling and crossover phenomena 

Let us now move from one-dimensional systems to quasi-one- dimensional systems by which we mean weakly 
coupled chains. A dimensionality crossover in low dimensional fermion systems is induced by a small coupling 
between chains. The nature of this crossover will depend on the kind of coupling involved. Weakly coupled chains 
correspond to so-called quasi-one-dimensional anisotropy, a situation realized in practice for electronic materials 
such as the organic conductors, 18, ^ lattice of spin chains, etc. In ordinary critical phenomena, scaling 



concepts were extended to anisotropic systems and have been successful in describing the new length scales 



FIGURE 1. One-fermion interchain hopping. ^ is the fermion coherence length along the chains. 



introduced by anisotropy; ^, |^ their use in the study of anisotropic fermion systems appears therefore quite 
natural. For example, deconfinement of quantum coherence for single particles can give rise to the restoration 
of a Fermi liquid component below some characteristic low energy scale. In addition, the deconfinement of pair 
correlations in more than one spatial dimension can lead to the emergence of true long-range order at finite 
temperature. |^ 

Let us consider first the case of interchain single-particle hopping (Fig. whose amplitude t± is small compared 
to the Fermi energy of isolated chains {tj_ <^ Ep)- Applying the extended scaling ansatz to the single fermion 
Green's function at | fc |= kp, and u!n=o — ttTjQ one can write 

G{T,k^) w DT~^XG{Bt^{k^)/T'^''). (1.19) 

i^(fcj^) = t±^ cosk± and k± is the transverse momentum of the Fermion spectrum (here the inter-chain distance 
d± has been set equal to unity), > is the crossover exponent which governs the magnification of the single- 
particle perturbation t±_ as ^ increases along the chains (Fig. |^). XG{y) is called a crossover scaling function 
and it determines transitory aspects of the crossover and the temperature range over which it is achieved. Here 
^g{^) = 1- When y <^ 1, the coherence of the fermionic system is confined along the chains and the physics 



is essentially dominated by the one-dimensional expression (1.7). Otherwise for y ^ 1, a different temperature 
dependence is expected to arise. The temperature scale for the single-particle crossover is thus determined by the 
condition y '--^ 1, that is 

T^-^y'^^ (1.20) 
Well below and on the full Fermi surface k — ]ip, one has 

G{T) ^ z{ti_)T~\ (1.21) 

where z{t\^ is the quasi-particle weight and 7 = 1 is the Fermi liquid exponent that restores the canonical 



dimension of the Green's function. In order to match (1.19) and (1.21) well below Tx^ one must have 



z{t^) = z^tl^^-''^'^-, (1.22) 

where z^o is a non universal constant. Applying scaling to the dynamical — coherent — part of the single-particle 
Green's function, the recovery of Fermi liquid behavior yields 

G'(k,w„)» . (1.23) 

with 7 = 1 again. 

Interchain coupling will affect the temperature dependence of pair correlations as well, which can undergo a 
dimensionality crossover and in turn an instability towards the formation of long-range order either in the Cooper 
or the Peierls channel at finite temperature. This possibility is of course bound to the number of components of 
the order parameter and the spatial dimensionality of the system. |^ The scaling behavior of pair correlations is 
therefore analogous to the one found in ordinary anisotropic critical phenomena. Properties of the system close 
to the true critical point at Tc can then be described by the usual reduced temperature 

T —T 

i--^- (1-24) 

The static pair susceptibility at some ordering wave vector Qo = (go, qi-o) will then diverge as t — *■ according to 

x{T)^C{g^)i-^, (1.25) 
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where 7 is the susceptibihty exponent. 

Applying the extended scahng hypothesis to the static susceptibihty in the superconducting or the density-wave 
channel, one can write 

x{T)r^ DT-'^X{Bg^/T'l'-^), (1.26) 

where X{y) is a crossover scaling function {X{Qi) = 1), is the intcr-chain coupling for pairs of fermions and 
0a;2 is a iwo-particle crossover exponent. Q Following the same arguments as in the single-particle case, the 
change in the exponents will occur when y 1, which allows us to introduce the crossover temperature for pair 
correlations 

^..^<?y^^^ (1.27) 

It is worth noting that the interchain coupling g±^ for pairs of particles differs from t^, and the two-particle 
crossover can occur above its single particle counterpart T^,, in which case the restoration of Fermi liquid behavior 
is absent. Instead only pair correlations deconfine, possibly leading to long-range order. 

If the scaling hypothesis holds, the expressions (1.25) and (1.14) have to match near Tc, which implies 

C{g^)^C^9l^'~'^'''^\ (1.28) 

where Coo is a non universal constant. Another prediction of the scaling hypothesis concerns the shift in Tc from 
zero temperatureQ 

Tc{9i_) ^ ^ 3^^==^ (1.29) 

When the system undergoes its two-particle dimensionality crossover, the temperature dependence of singularities 
near Tc is governed by the correlation length 

(1.30) 

where is the coherence length at T ^ T^;, which fixes a short anisotropic length scales for the onset of critical 
fluctuations described by the exponent v. At Tc, one has the critical form 



which according to (1-25) leads to the Fisher relation 

7=(2-77)i>. (1.32) 



3 Free fermion limit 

Before facing the task of calculating the exponents and various scaling quantities of interacting fermion models, 
it is instructive to first analyze the properties of the free fermion gas model. We will consider first the purely ID 
system and then introduces interchain single particle hopping and see how they actually satisfy the above scaling 
relations. 



3.1 One dimension 

Let us consider a one-dimensional system of free fermions coupled to source pair fields h for which the partition 
function is given by a functional integral over the anti-commuting (Grassman) fields ^p (see references for 
reviews on Grassman fields) 

Zo[h*,h] = J J Vi^*Vtlj e'^oir M+Snir ,i',h-\h] 

- exp{^[G°(:fc)]-V;,a(^)^P.a(fc) + E[OM®^M(9) + c.c]}. (1.33) 

k.a li,q 
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Here 5*0 is the action of the system in the absence of fields and 

G°(fc) = ~{^(k)r(k))o 

[iujn - (1.34) 

is the free fermion propagator in the Fourier-Matsubara space where k — {k,ujn). The fermion spectrum is 
linearized with respect to the Fermi points ifej?, namely 

ep{k) — vp{pk — kp), (1.35) 

where p — ± refers to right and left (— ) going fermions and vp is the Fermi velocity. A natural cutoff on 
the bandwidth can be imposed on this continuum model by restricting the summations on the wave vector k of 
branch p in the interval pkp + ko > k > pkp — kg. The cutoff fco ^ where a is the lattice constant on the 
band wave vector, then leads to a finite bandwidth Eq = 2vpko = 2Ep which we will take equal to twice the 
Fermi energy Ep = vpkp. 

In one dimension, the fermion gas is unstable to the formation of pair correlations in both Cooper and Peierls 
channels. This shows up as infrared singularities for pair susceptibilities, which can be calculated by adding Sh to 
the action, which linearly couples source fields to fermion pairs. In the Peierls channel (/x = fip), the particle-hole 
fields 

O^p=o® = E ^-.-(^ - Q)^+Ak) (1.36) 

and 

0,,M^ = E ^-,"(^ - q)<^f,^+Ak) (1.37) 

correspond to charge-density- wave (CDW: — 0) and spin-density- wave (SDW:^p ^ 0) correlations respectively 
near the wave vector 2kp, where q = {2kp + q,u!m) , o'pj,=i.2,3 are the Pauli matrices and L is the length of the 
system. In the Cooper channel (/i = /ic)i the particle-particle fields 

k,a 

and 

O^^Mq) = E V'-,-a(-fc + q}af^P+Ak), (1.39) 

correspond to singlet (SS: = 0) and triplet (TS: /ic 7^ 0,) superconducting channels respectively, where 
q= {q,LUm)- The susceptibilities in zero /i^ in both channels are defined by 



,^ 1 S^Zo[h*,h] 



Z„ SKiq)Sh,{q) 



h=0 



- -{O^iq)0;{q})o, (1.40) 
and obtained from the calculation of the statistical averages 

X,{q) - - Z^J I VrVi^ 0,{q}0;{q) e''"^^' '^l (1.41) 

in the Peierls and Cooper channels. Thus in the Peierls channel, the dynamic susceptibility (Figure ^a) close to 
2kp is given by the standard result 

Xt^Ao) = 2'^J2'^-(k~2kp^q,LOn~LUrn)G'Ak,LUn) 
k 



D 



D 



D 



FIGURE 3. Energy scales for the ID fermion problem. Between the Fermi (Ep) and thermal energy (T), there is no 
characteristic energy scale. 



This logarithmic divergence can also be seen as a power law with a vanishing exponent, that is 



X^^iT)^ hm -7V(0h-i 
7— >o+ 



T \ -7 



I.IS^F 



) 



(1.46) 



According to the relation (1.18), this should imply 77 —s- 2, namely a logarithmic singularity in q. This is indeed 
consistent with the asymptotic behavior of Xi-iil) at T = and ujm = 0, which is found to be 



\VFq\ 
2Ef 



(1.47) 



in agreement with ( 1.15 ). 

The origin of scale invariance can be further sharpened if one looks at the decay of the fermion coherence at 
large distance. Consider the Green's function 



T 



^ikx—iuJnT 
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27ri smh[{x + ipvpT)/^] 

where 



(1,48) 



f = ^ (1.49) 

is a characteristic length scale for the quantum coherence of each fermion and in fact the de Broglie quantum 
length. At low temperature for a <C x ■f, one finds the power law decay of coherence as a function of distance 

G°(x) ~ 1/x, (1.50) 

which signals homogeneity or scaling. Reverting to ( |l.6| ), we verify that = in the free fermion limit. For a; 3> f , 

Glix) - e-"/« (1.51) 

there is no homogeneity and the coherence is exponentially damped with distance. The absence of energy scale 
between Ep and T thus corresponds in real space to the absence of length scale between ^ and a. Rewriting the 
elementary the Cooper and Peierls susceptibilities as 

Xp(T)«-7V(0)ln^, (1.52) 
a 

one then sees that the quantum spatial coherence of particles becomes an essential ingredient in building up 
logarithmic singularities in pair correlations of the Cooper and Peierls channels. This can be directly verified by 
the decay of pair correlations in space and time. From the definition (1.13) and the use of Wick theorem for free 
fermions, one obtains 



. 1 1 



27r2^2 sinh[(a; + ivpr)/^] sinh[(a; — ivpr)/^] ' 
where the upper (lower) sign corresponds to the Cooper (Peierls) case. For a ^ x ^ one finds 



(1.53) 



e 



iqox 



Xt^ix) oc (1.54) 



in agreement with ( 1.14 ), in which d — 2 and rj — 2. At large distance, when a; S> ^, one has 
where the effective coherence length for pairs is half the one of a single fermion. 



3.2 Interchain coupling 

A very simple illustration of a dimensionality crossover for the single-particle coherence can be readily given in 
the free fermion case by adding an interchain hopping term to the the free action, which becomes 

k,p,a (^.J) 

= EE[^p('^'^")]"V;a(k,u;„)Vp,a(k,c.„) (1.56) 

k,uj„ p,a 

where t±^ is the single-particle hopping integral between nearest-neighbor chains i and j. Considering a linear 
array of iVj^ chains, the free propagator that parameterizes Sq in Fourier is 

GO(k,c^„) = (iL^„-£;p(k))"\ (1.57) 

where 

Ep{\i) = ep(fc) - 2t^ cos fc_L (1.58) 
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is the full fcrmion spectrum. The propagator can be trivially written in the extended scaling scaling form (1.19) 
with 7=1, where the scaling function takes the simple form 

XG{v) = {l + vr\ (1.59) 

with y = Bt± cos{k±)/T and B = {iir /2)^^. This implies (jj^i = 1 for the single-particle crossover exponent. This 
value of is found in the anisotropic free field theory of critical phenomena. Q The Fourier transform leads to 
single-particle fermion correlation function 

G°(a;,ni,r) =G°(a;,r)e-'""^J„^(a;/C.), (1.60) 

where n±_ is the transverse separation expressed in number of chains. The crossover scaling function is then 
embodied in the Bessel function J^^ ix/S,x) of order nj_, which takes sizable values when the intrachain distance x 
reaches values of the order of the characteristic length scale — VF/{2t±) needed for coherent hopping from one 
chain to another (Figure p. For x ~ ^, the crossover condition is ^ ^ ^x, and this leads once again to T^i ~ t±, 
namely (jj^^ — ^■ 



4 The Kadanoff- Wilson renormalization group 

4-1 One- dimensional case 

When the motion of fermions is entirely confined to one spatial dimension, particles cannot avoid each other 
and hence the influence of interactions is particularly important. The absence of length scale up to ^ found in 
the free fermion gas problem for both single and pair correlations will carry over in the presence of interactions, 
which couple different energy scales. The low energy behavior of the system in which we are interested will be 
influenced by all the higher energy states up to the highest energy at the band edge. Renormalization group 
ideas can be used to tackle the problem of the cascade of energy or length scales that characterizes interacting 
fermions. This will be done in the framework of the fermion gas model in which the direct interaction between 
fermions is parametrized by a small set of distinct scattering processes close to the Fermi level (Figure ^) . |^ This 
is justified given the particular importance shown by the Cooper and Peierls infrared singularities for electronic 
states near the Fermi level. We will consider the backward (gi) and forward (32) scattering processes for which 
particles on opposite sides of the Fermi level exchange momentum near 2kF and zero respectively. When the band 
is half-filled, Akp coincides with a reciprocal lattice vector, making possible Umklapp processes, denoted by 33, in 
which two particles can be transferred from one side of the Fermi level to the other. Thus, adding the interaction 
part Si to the action, the partition function in the absence of source fields takes the form 



JJvrv^P exp{^[G«(fc)]-V;,„(fc)V'p,a(fc) 



k.a 



+ nVpT/L V+.aiC^l +5)V'l,a2(fc2 - g)V'-,a3(fc2)V'+,a4(fcl) 

{p,k,q,a} 

- nVFT/2L ^ 33 V'p,a(^l +P9)V'lp,a'(fcl -M+pG)V'-p,a'(A2)-(/'p,a(^l)| (1.61) 



where 



9{a} — 9l^aia^^a2a4 .92^aiQ4^Q2CK3 (-^•^^) 

and G = {4:kp, 0) is the reciprocal lattice vector. The parameter space of the action is then 

t^S = {Gp, 91,92,93)- (1-63) 



where all the couplings are expressed in units of t:vf- 



FIGURE 4. Backward (ffi), forward (32) and umklapp (gs) couplings of the 1-D fermion gas model and the corresponding 
diagrams. The open (full) circle corresponds to the generic vertex part for backward and forward (umklapp) scatterings. 



The idea behind the Kadanoff- Wilson renormalization group is the transformation or renormalization of the 
action S following successive partial integration of -tp's of momentum located in the energy shell ^Eo{i)di on both 
sides of the Fermi level, where Eo{£) = Eoe~^ is the effective bandwidth at step £ of integration and d£ 1. The 
transformation of S from £ to i + d£, that keeps the partition function invariant is usually written as 

e^(^+'^^) I?V*25V'e^'^*'^''+"^ (1-64) 
where A{£) corresponds to the free energy density at the step £. The integration measure in the outer shell is 

vrv^p= Yi drpjk)d^PpAk)- (1-65) 

'p,a,{k}' 

and {k}' — {fc}'{(jj„}. Here {k}' corresponds to the momentum outer shells 

fcoe"^ + kp > k > koe^^^"^^ + kp {-kQC^^^'^^ - kp > k > -k^e^^ - kp) 

above the Fermi level and 

-fcoe~^~''^ + kp > k > -fcoe"^ + kp (-fcoe"^ - kp > k > -fcoe~^~'"' - kp) 

below for right (left) moving fermions, while {w„} covers all the Matsubara frequencies. The remaining inner shell 
(<) fermion degrees of freedom are kept fixed. The parameters of the action denoted as fis transform according 
to 

RMt^s{e) = t^sie + d£). (1.66) 

where 

fisii) = {z{£)G% zi(£)gi, Z2(%2, Z3(%3) (1.67) 

contains the factors z and 2^=1,2.3 that renormalize the parameter space of the action at ^.[Q The flow of 
parameters is conducted down to £t = InEp/T corresponding to the highest value of £ at which scaling applies 
as discussed in § II. This leads in turn to the temperature dependence of the parameter space /is(T'). It should 
be mentioned that the Kadanoff transformation will generate new terms which were not present in the action at 
^ = 0. In principle, the relevance of these can be seen explicitly from a rescaling of energy and field 'density' after 
each partial integration. Thus performing an energy change by the factor s = e'^^ > 1, one has 

a;„ = sujn- (1.68) 
The fermion fields in the action will then transform according to 

= z-5s-^(*), (1.69) 

which in its turn implies that the rescaling of interaction parameters is given by 

gl^g,z,z^, for i-1,2,3. (1.70) 



14 C. Bourbonnais, B. Guay and R. Wortis 



This contains only anomalous renormalization factors indicating that the interaction parameters are all marginal. 
However, if interactions between three, four or more particles were added to the bare action, the amplitude 
of these would decay as 1/s to some positive power, and would therefore be irrelevant in the RG sense. A 
different conclusion would be reached, however, if these many-body interactions were generated in the flow of 
renormalization. In this case, the amplitudes of interactions becomes scale dependent and such couplings can 



remain marginal or even become relevant (see § 4.3 and § 5.1). Besides this use, rescaling for our fermion problem 



is not an essential step of the RG transformation. We will consider the renormalization of us with respect to £ as 
a transformation that describes an effective system with a reduced band width Eq {£) or a magnified length scale 
vp/Eoii). 



In practice, the partial integration (1.64) is most easily performed with the aid of diagrams. We first decompose 



the k sums in the action into outer and inner shells momentum variables 

E-E+E- (1-71) 

{k} {k}' {fc}< 

This allows us to write 

Sir ,^p] = s[r,^]< + s[r, ^, r (1.72) 

where S'['0*,'0]< of the action with all the ip's in the inner shell, whereas 

4 

s[r,^,r,^]^So[r,^P] + Y,siAr,^,r,i^] (1-73) 

consists in a free part in the outer momentum shell and an interacting part as a sum of terms, Sj^i, having 



i = 1 , ... 4 ip's in the outer momentum shell. Their integration following the partial trace operation ( 1.64 ) is 
performed perturbatively with respect to So[ip*,4']- Making use of the linked cluster theorem, the outer shell 
integration becomes 

4 

J vrv^e^^ [s[r,^]< + E V"! (( E siA^,^. r,^] )">o,c} (1-74) 



where 

(.•.•)o,c = Z^' J I' Vrv^ ( ) e^-'I'^*''?! (1.75) 

is a free fermion average corresponding to a connected diagram evaluated in the outer momentum shell and 

Zo^ JJ I?V;*2?Vie^°['^*''^l, (1.76) 
is the outer shell contribution to the free partition function. 



4-2 One-loop results 

The KW method may be used to obtain low order results presented in previous review papers, |l^ and no 
particular technical difficulties arise. The following two subsections present material which has already been 
developed using the parquet summation]^, ^ and the first order multiplicative renormalization group. ||lO|, ^s] , pT[ 
Nevertheless, it is useful to supply the basic properties of the logarithmic theory of the ID fermion gas model 
using this method. 



D 



FIGURE 5. Recursion relations for backward (gi) and forward (32) scattering at the one-loop level. The slashed lines refer 
to a particle or a hole in the outer shell. 

With an obvious rearrangement, the flow equations may be written 

dgi 2 

^(2g2-5i)=0. (1.81) 
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The solution for gi {£) is found at once to be 



= (1-82) 



showing that backward scattering is marginally irrelevant (relevant) if the bare coupling gi is repulsive (attractive). 
As for the combination 2g2 — gi, it remains marginal, an invariance that reflects the conservation of particles on 
each branch. ^ Another feature of these equations is the fact that the flow of gi is entirely uncoupled from 
the one of 2g2 — gi- This feature, which follows from the interference between the Cooper and the Peierls channels, 
gives rise to an important property of ID interacting fermion systems, namely the separation of long wavelength 
spin and charge degrees of freedom. This is rendered manifest by rewriting the interacting part of the action as 
follows!^ 

Si[ip*,ip]e = -TTVF{2g2- gi)^Pp{q)P-p{-q) + nvp gi{e)^Sp{q) ■ S^p{-q) (1.83) 

p. 9 



where the long- wave length particle-density and spin-density fields of branch p are defined by 

Pp® = ^ ^p,a(fc + 5)V'p,a(fc) 

Sp(5) - E rpAk + q)<f''^i^pAk). (1.84) 

a,{/e}< 

The spin-charge separation is preserved at higher order and is a key property of a Luttinger liquid in one 
dimension. In the repulsive sector 171 > 0, the spin coupling goes to zero, while the forward scattering 
52 ~^ ~^ 92 — 5i/2 flows to a non universal value on the Tomanaga-Luttinger line of fixed points gl — 
(Figure §). 

For an attractive interaction 171 < 0, an instability occurs in the the spin sector at i„ defined by giia = ^1- 
Identifying £ with £t = InEp/T, the temperature dependence of the Peierls or Cooper loop, the characteristic 
temperature scale for strong coupling is 

T, = EFe-^/^S'^ = ^A„, (1.85) 



where A„ is defined as the spin gap. According to (1.83), backward scattering can be equated with an exchange 



interaction between right- and left-going spin densities which is antiferromagnetic in character when gi is at- 
tractive. Therefore strong attractive coupling signals the presence of a spin gap A^. Although the divergence in 



(1.82) is an artefact of the one-loop approximation, the flow of gi{£) to strong coupling is nevertheless preserved 
at higher order. This is shown by the fact that the system inevitably crosses the so-called Luther-Emery line at 
5i(^le) = —6/5 (Figure ||), where an exact solution using the technique of bosonization confirms the existence 
of a spin gap.Q |o) 

Umklapp scattering at half- filling 

At half-filling, there is one fermion per site so that Akp = 27r/a coincides with a reciprocal lattice vector G enabling 
Umklapp scattering processes denoted by 53 in Figure |[ The presence of g^ processes reinforces the importance 
of the 2fci? density-wave channel compared to the superconducting one. Flows will be modified accordingly. Thus 
at the one-loop level, the outer shell corrections coming from ^{{81^2)^)0.0 will include additional contributions 

5/^2 -^=» S'/,2|g^^^o + (■0+V+i/'-'0- + perm. -f c.c) 

+ {ij*+^l^-^- + ^*+^*+^-^- + c.c) (1.86) 

involving outer shell fields in the Peierls channel and pairs of particles or holes on the same branch. At the one- 
loop level only the former gives rise to logarithmic corrections, whereas the latter will be involved in higher order 
corrections (§ 4.3). The former contributions, along with ( 1.78| ), will lead to the recursion formulas 54=1,2.3 (^+c?^) = 



Zi=i,2,3{d£)gi=i.2.3{^), which can be written in form 

gi{£ + d£) = gi{£)- 2gl{£)Ip{d£), 



FIGURE 6. Flow of coupling constants in the gig2 plane at the one-loop level. In the repulsive region gi > 0, the system 
has a line of fixed points on the Tomonaga-Luttinger line — 0. In the attractive sector gi < 0, the fiows is to strong 
coupling. The dashed line corresponds to the Luther-Emery line at gi — —6/5. 



(2g2-gi){£ + M) = {2g2 - gi){e) + 2gl{i)Ip{M), 

gsie + di) = g3{£) + 2{2g2-gi){i)g3{i)Ip{di). (1.87) 

The flow of the coupling constants are then governed by 

dgi _ _ 2 

di ~ 
d 2 

^(2g2-.gi) = 93, 

^ = 33(252-51), (1-88) 

which corresponds to the diagrammatic parquet summation of Dzyaloshinskii and Larkin. |^ By combining 
the last two equations, the flow is found to follow the hyperbolas described by the renormalization invariant 
C = [{2g2 — 5i)(^)]^ — ffK^)- Umklapp scattering is then entirely uncoupled from 51 and will only affect the charge 
sector. When 51 — 252 >| 53 |, 53(^) scales to zero and is marginally irrelevant, while 252 — gl is non universal. 
The excitation of charge degrees of freedom are then gapless in this case. The attractive Hubbard model in which 
51,2,3 = U < falls in this category. 

Otherwise, the flows of 53 and 252 — 51 evolve to strong coupling and both couplings are marginally relevant. 
There is a singularity at £p corresponding to the temperature scale 

Tp^EFe-'/V\^\ =^Ap. (1.89) 

Following the example of the incommensurate case, the couplings inevitably cross the Luther-Emery line at 
5i — 252 = —6/5, at which point the problem can be solved exactly and the existence of a charge gap Ap is 
confirmed, The special case of the repulsive Hubbard limit at half-filling where 51,2.3 = U > falls on the 
separatrix 252 — 51 =| 53 | (Figure 0). In this case the solution of (E8S) leads to a singularity at Tp = £'i?e^^/l3-'l . 



A charge gap Ap is consistent with the characteristics of a ID Mott insulator found in the Lieb-Wu exact solution. 



4-3 Two-loop results 

At the one-loop level, corrections to the action parameters fis resulting from the partial trace are confined to 
the scattering amplitudes 5i=i,2.3(^) and aside from a renormalization of the chemical potential, the flow leaves 
the single particle properties unchanged. At the two-loop level, outer shell logarithmic contributions will affect 
both the one-particle self-energy and the four-point vertices. This yields an essential feature of the ID electron 
gas model which is the absence of quasi-particle states. 

One-particle self- energy 

Let us consider first the one-particle part. It is easy to see that the term ^ ((fS*/, 3)^)5 at n = 2 has the form of 
a one-particle self-energy correction. However, if we strictly follow the classical Wilson scheme, the constraint of 



FIGURE 7. The flow of couplings in the charge sector. The condition gi — 2g2 > \ ga \ corresponds to weak couphng where 
Umklapp scattering is a marginally irrelevant coupling. Otherwise, both gi — 2g2 and g^ axe marginally relevant. The 
Luther- Emery line is located at gi — 2g2 = —6/5. 



having all internal lines in the outer shell leads to a vanishing contribution. It becomes clear then that one-particle 
corrections can only be logarithmic at the two-loop level if the internal lines refer to different momentum shells. 
In the classical Wilson scheme, this possibility actually results from a cascade of contractions on the three-particle 
interaction (Figure ||). Such interactions, although absent in the initial action, are generated along the RG flow 
and enlarge the parameter space. 



In order to see how these interactions are generated from the partial trace (1.74) and under what conditions 
they give rise to logarithmic corrections at the two-loop level, it is useful to label the steps of the outer shell 
integration by the index j, for which i?o(^j") = EQe~^^~'^'^'^^ is the scaled bandwidth at j > 3. Thus at each step 
j — 2, the contraction ^{Sj ^)q ,, yields a three-particle term in the action (Figure |^a), which is denoted SS^^\ Its 
explicit expression for the forward and backward scattering is given in Eq. ( p77| ) of § The analo gous expression 
^^\u Umklapp with A J as the scattering amplitudes is given in (|182|) of § Owing to the presence of the 



\u 

internal line, the rescaling of energy and fields shows that the amplitudes and Aj are marginal, At the 
next step of partial integration j — 1, two fermion fields of these three-particle terms can in turn be contracted 
in {SSi^\)Q ,,, from which we will retain an effective two-particle interaction with two internal lines tied to the 



adjacent shells j — 2 and j — 1 (Figure H-b and (178)). These contributions are non- logarithmic and differ from 



those already retained in the one-loop calculation in §4.2 for which both internal lines were put in the same outer 
shell. The repetition of these contractions is carried out by a sum on j up to — 1 (Figure |^c) , which yields 

^S'x = J2f~\SSi^,l)o,c i^S^^ = J2f~\^si^!u,2)o,c for Umklapp) as shown in (|l7|) ((|l||) for Umklapp). At the 
next step, A^, of partial integration, two more fields arc in turn contracted which lead to a sum of one-particle 
self-energy corrections for a k state at the top of the inner shell. [5^ Since each term of the sum refers to couplings 
at different steps j, the correction is non-local. 

In the following, we will adopt a local scheme of approximation, which consists in neglecting the slow (log- 



arithmic) variation of the coupling constants in the intermediate integral over momentum transfer q in (181) 
and evaluate them at £ (Figure ^d). It is worth noting that the local approximation becomes exact for cou- 
plings or combination of couplings that are invariant with respect to £. This turns out to be the case for the 
Tomanaga-Luttinger model when only g2 coupling is present (gi = = 0), or for the combination 2g2 — gi when 
33 = 0. Neglecting small scaling deviations, the end result given by (181) (( |l84| ) for Umklapp) is logarithmic and 



contributes to a 1/z renormalization factor of [Gp{k,ujn)] ^ for k located at the top of the inner shell. The flow 
equation for z{£) is 

^ Inz = -^[(232 - 9i)' + igl + 233^]. (1.90) 
The solution as a function of £ can be written as a product z — z„Zp of renormalization factors, in which 

z,(£)=exp{-^^ (1.91) 

for the spin part and 



Zp{l) = exp { - 1 ^ [(252 - g^f + 2gl]d£'] 



(1.92) 



FIGURE 8. The cascade of contractions contributing to logarithmic one-particle self-energy corrections in the local ap- 
proximation for the RG flow: (a) backward and forward scattering contributions to the three-particle interactions; A"*" (left) 
and A~ (right) from |(S|,i)o,c a^t the step (j — 2) where j > 3; (b) second contraction in the outer shell at the step j — 1 
with two internal lines on successive energy shells (only shown for A^) ; (c) the summation of diagrams given in (b) up to 
j = N — 1 ; (d) outer shell one-particle self-energy correction at the step A'^ [Ndi — e^^"^^) for the top inner shell state k of 
the +kF branch in the local approximation. Here a crossed (slashed) line refers to fermions in the interval [0,^] dl]); 
(e) the result of a similar cascade of contractions for Umklapp scattering. 



for the charge. Prior to analyzing the impact of these equations on one-particle properties, to work consistently 
at two-loop order, we must proceed to the evaluation of the two-loop contribution to the coupling constant flow. 



Four-point vertices 

To obtain the two-loop corrections to four-point vertices within a classical KW scheme, we come up against 
difficulties similar to those found for the one-particle self-energy. We can indeed verify that contractions of the 
form {S'j 2)0,0 8.t n — 3, which would lead to two-loop corrections for the four-point vertices, are actually 
vanishingly small if all internal lines are in the same momentum shell. Here again non-zero contributions can 
only be found if the internal lines refer to different shells, which implies that the contractions must be made 
starting from interactions that involve more than two particles (Figure |9|).[p^ Here we skip the technical details 
of the calculation and only sketch the main steps, which run parallel to those of the one-particle self-energy. The 
cascade of contractions starts from the four-particle interaction resulting from 2'S'/,i5'/,i)o c n = 3, which 
involves the contribution of the Landau channel. The corresponding diagrams are shown in Figure ^a at 
the step j — 2. For the required contractions to be possible, the momentum transfer must be zero. At the next 
step, two more fields are contracted and a line is put in the adjacent shell j — 1 (Figure ||-b). As a function of j, 
these corrections add up as shown in Figure ^c to lead to effective three-particle interactions at the step N — 1. 
At the last step, N, of the cascade, one arrives at a two-loop correction for the four-point vertices (Figure ^d). 
Following the example of one-particle self-energy, these corrections are non-local in the coupling constants. In 
the local approximation scheme, one finds an outer shell contribution that is 0{g^d£) (Figure ^d). The two-loop 
contributions involving Umklapp scattering are shown in Figure 0-e. Finally, the outer shell contributions to the 



FIGURE 9. The cascade of contractions contributing to the four-point vertices in the local approximation of the RG flow: 
(a) backward and forward scattering contributions to the generation of four-particle interactions from | (S'f'25/,iS/,i)o,c 
at the step j ~2 {j > 3) ; (b) second contraction in the outer shell at the step j — 1 ; (c) the summation of diagrams given 
in (b) up to j' = A'' — 1 as effective three-particle interactions; (d) two-loop outer shell correction to four-point vertices in 
the local approximation at the step Nd£ 1 + dl. Here a crossed (dashed) line refers to fermion lines in the interval [0, P\ 
+ dl])\ (e) the result of a similar cascade of contractions for Umklapp scattering. 



Now, combining these with the one- loop results of ( 1.88 ) and performing the rescahng of the fields, one has 

^-l/2^(*) ^ ^(*) 

^is = {Gp,ziz'^gi,Z2z'^g2,Z3z'^g3, . . .), (1.94) 

which renders explicit the rcnormalization of the couplings gi. Neglecting scaling deviations at small £, the 
corresponding scaling equations are 

dgi _ _ 2 _ 1 3 

-(252-51) = 53[l-^(2<?2-.gi)], 

^ = 53(252 -5i)[l-J(252 - 51)1-^53', (1-95) 

which reproduces the old results of the multiplicative renormalization group method, [jll], Q We see that the spin 
coupling 51 is still uncoupled to the charge couplings 252 — 51 and 53. Furthermore, the one- loop singularities. 
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previously found at T^^^^p, are removed and replaced by strong coupling fixed points for the couplings as ^ ^ cxi. 
For the interesting cases already analyzed at the one loop level, we have first the attractive one where gi < and 
91 ~ 2(72 > \ 93 I J for which gl —2, (73 0, while gl — 2g2 is strongly attractive but non universal; second, in the 
repulsive case when gi > and gi — 2g2 < | 33 |, for which Umklapp is marginally relevant with 53 — > 2, 1 
and gl — s- 0. Since in both cases, the couplings cross the Luther-Emery line, the scale for the gap remains of 
the order of the one-loop temperature scale T^, in both sectors. 



The integration of (1.90) at large £^ ii, then leads to the power law decay 

z{e)^{Eo{e)/A,f. (1.96) 

The exponent evaluated at the fixed point in the presence of a single gap can be separated into charge {6p) and 
spin (6*0-) components: 

r ^e;i2g;~gig;) + e;igi) 

~ 3/4. (1.97) 

The case when there is no Umklapp and gi > corresponds to the Tomonaga-Luttinger model at large £. The 
couplings remain weak and one finds the power law decay 

z{£)^D„{i)(^Eoie)/Eoy\ (1.98) 

with 

^i2g2-gi)^ (1.99) 

and ^ 

D,{i) = exp ( - ^ ^ 9!{ndl') ■ (1.100) 

Here the exponent 9p is non-universal. |ll] A contribution from the spin degrees of freedom appears through the 
weak transient Da{tj- The physics of this gapless case is that of the Luttinger liquid. 

The factor z is the quasi-particle weight at Fermi level. In all cases, it vanishes at ^ 00 or T ^ indicating 
the absence of a Fermi liquid in one-dimension. Actually the same factor also coincides with the single-particle 
density of states at the Fermi level 

N*{£) = N{Q)z{t), (1.101) 

which vanishes in the same conditions. 

The connection between the one-particle Green's function computed at ^ = and that computed at non-zero 
€ can be written in the form ^ ^ 

Gp(fc, tis) = z{i)Gp{k, fisii)). (1.102) 

When £ ^ 00, the number of degrees of freedom to be integrated out goes to zero, thus Gp{k, Hs{£)) zG^{k 
kp, LUn) and as a function of T — )■ 0, one finds 

Gp{kF,L,r.=0,^^s)~CT~\ (1.103) 



where, consistent with the scaling ansatz (1.7), 7 = 1 — 6* and C — \/{i'KEp). If we identify the vanishing energy 
^i?o(^) in z{t) with vp{pk — kp) close to kp, we have at T = 

Gp{k,^is)oi\pk-kp\-^+\ (1.104) 

which agrees with the scaling expressions ( |l.5| - |l.6[) showing the absence of a Fermi liquid. 

Remarks 

At this point, casting a glance back at the RG approach described above, a few remarks can be made in connection 
with the alternative formulation presented in Refs. ||l^, ^ In the classical Wilsonian scheme presented here, a 
momentum cutoff fco (or bandwidth Eq) is imposed on the spectrum at the start and all contractions done at a 
given step of the RG refer to fermion states of the outer shell tied to that particular step. In this way, we have 
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seen that the two-loop level calculations require a cascade of outer shell contractions starting from many-particle 
interactions that were not present in the bare action. The cascade therefore causes non-locality of the flow which 
depends on the values of coupling constants obtained at previous steps. Thus obtaining logarithmic scaling at 
higher than one-loop forces us to resort to a local approximation. 

In the formulation of Refs. , the Kadanoff transformation at high order is performed differently. In effect, 

a cutoff fco is not imposed on the range of momentum for the spectrum in the free part of the action 5*0 but it 
only bounds the A:-summations of the in teraction term Si . The latter summations, of which there two in both of 
the interaction terms (as in Eqn. ( 1.6l| ) ), are involved in the ultraviolet regularization of all singular diagrams, 



while the remaining summation on the momentum transfer q is treated independently. In this way, the interaction 
term Si jij), ■0*, ■0] with fields to be integrated out in the Kadanoff transformation is defined differently than in 



(1.73). It has no Si^i term with a single field in the outer shell, while Si_3 and Si^i contain a momentum transfer 



summation corresponding to the high-energy interval | vpq |> Eo{£)/2. The generation of new many-particle 
terms via Si,i is therefore absent, and the evaluation of 5((»5/^3)^)q ^ for the one-particle self-energy correction, 
directly yields (aside from some differences in the range of momentum transfer integration and in the scaling 



deviations) the end result of the local approximation of § 4.3 at the two-loop level. A similar conclusion holds for 
the four-point vertices when the calculation of ^{S] ^Sf 2)0.0 is performed. Despite its simphcity, however, this 
alternative formulation is less transparent as to the origin of logarithmic scaling in the ID fermion gas model at 
high order. Moreover, the presence of a summation over a finite interval of momentum transfer for each Kadanoff 
transformation removes from the fermion spectrum its natural cutoff fco and stretches the standard rules of the 
Wilson procedure. 

4.4 Response Functions 

We now calculate the response of the system to the formation of pair correlations in the Peierls and Cooper 
channels. In the Peierls channel, density-wave correlations with charge or spin buildup centered either on sites or 
bonds differ at half- filling due to the presence of Umklapp scattering. Actually, because 4fci? is a reciprocal lattice 
vector, density waves at +2kp and —2kF turn out to be equivalent. In a continuum model, one can thus define 
the following symmetric and antisymmetric combinations of pair fields 

0^±{q)^^{0;^{q)±0^,{q)), (1.105) 

corresponding to site (+) and bond (— ) (standing) density-wave correlations. 

We consider the linear coupling of these, as well as those previously introduced for the Cooper channel in 



(1. 38-1. 39), to source fields /i^; the interacting part of the action Si then becomes 

Si[r,i^,h*,h] = si[r,^] + Sh[r,i^], (i.ioe) 

where 

=E('^M(5)0;(g) + c.c). (1.107) 



Following the substitution of the above Si in the action (1.72), the KW transformation will generate corrections 
that are linear in hp, and also terms with higher powi 
terms are kept so that the action at £ takes the form 



that are linear in hp, and also terms with higher powers. In the linear response theory, only linear and quadratic 



S[r,i',h*,h]e = S[r,i^]i + Y.{4(^Ki^Op{q) + C.C. 



where z^{£) is the renormahzation factor for the pair vertex part in the channel /x (Fig. p^a). The term which is 
quadratic in the source fields corresponds to the response function (Fig. |lO|-b) 

Xp{£) = -{nvp)-' f x^{£')d£\ (1.109) 
"'0 
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FIGURE 10. Diagrams of the flow equations of the pair vertex part (a) and susceptibility X/j (b) in the channel fi. 

In the gapless regime for both spin and charge when gi > and gi — 2g2 >\ gs |, the ^ =TS and fi —SS 
responses are singular 

X.iT) -{§^y''. (1.115) 
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FIGURE 11. Phase diagram of the ID interacting Fermi gas in the presence of Umklapp scattering. 



To do so, it is useful to restore the wave vector dependence of Xn by evaluating the upper bound of the loop 
integral in ( 1.109| ) at 



^-.%,T) = lnHi^+ VXl/2)-i 



V^(l/2 



+ i 



47rr 



c.c 



(1.116) 



which coincides with the complete analytical structure of the free Peierls and Cooper loop in ( |l.42 ) and (1.44) in 
the static limit, respectively. At low temperature when | vpq |S> 47rT — > 0, we can use (1.47) and obtain 



Xt^iq) oc q or Xt^i^) cc x^" ^, 

which are consistent with the scaling relation 2 — 77^ = 7^ at = 1 (Eqn. ( 1.1^ )). 



(1.117) 



5 Interchain coupling: one-particle hopping 



Let us now consider the effect of a small interchain hopping (Figure |l|) in the framework of the KW renormalization 
group. The bare propagator of the action at ^ = is then replaced by Gp(k, ct;„) given in (1.57). The parameter 
space for the action at ^ = then becomes 

MS = (G'p(fc,w„),tj_, 51, 52,53)- 



Now applying the rescaling transformations ( 1.68 - L69| ) in the absence of interaction, the rescaling of ID energy 
leads to the transformation 

t'j_ = stj_, (1.118) 

which means that the canonical dimension of t± is unity and it is therefore a relevant perturbation. An estimate 
of the temperature at which a dimensionality crossover in the coherence of single-particle motion occurs can 
be obtained by replacing s with the ratio ^/a and setting t'^ ~ Ep. We find T.j;i ^ tj_ in agreement with the 
results of § 3.2. How this is modified in the presence of interactions will depend on the anomalous dimension 9 
which the fermion field acquires under ID renormalization. When the partial trace is performed at low values of 
i, one-particle self-energy corrections are essentially ID in character and are governed by ( 1.90| ). At step £, the 
effective 2D propagator then takes the form 



Glik,uj„) = z{i)[Gl-\k) + 2t^z{i)coskj 



(1.119) 



FIGURE 12. Crossover temperature variation as a function of the one-particle self-energy exponent 8. 



Therefore the power law decay of the quasi-particle weight and the density of states (cf. Eqn. ( |1.10l| )) along the 
chains leads to a renormalizcd hopping amplitude ztj_ which decreases as a function of i. The condition for the 
single particle crossover temperature T^i is then renormalizcd, becoming z{T)tj_^/a ~ Ep, which implies 



■Ef' 



where = 1 — 6' is the single particle crossover exponent. ||55 

behaves as a Fermi liquid. The effective bare propagator (1 _ , „ , . „ 

has the Fermi liquid form (1^) with a weight that is consistent with the extended scaling expression (|l.22|) with 



(1.120) 

In this simple picture, well below T^i , the system 
for Ep[\<i) < ^Eo{£xi), where 4i = InEp/Tx 



7 



1 



7 



1 and the non universal constant 



E' 



As long as 6* < 1, i_L remains a relevant variable and T^i is finite. For 9 = 1 and 6* > 1, however, it becomes 
marginal and irrelevant respectively. In these latter cases, transverse band motion has no chance to develop and 
fermions are confined along the chains at all temperatures. In the KW framework, conditions for ID single-particle 
confinement are likely to be realized in the presence of a gap since in this case we have seen that 9^1. |5^, ^ 
Single-particle confinement is also present in the absence of any gaps, as in the Tomanaga-Luttinger model, where 
91 — 93 = a-iid 9^1 when | 32 | is close to unity js^, ^ and other interacting fermion models |59, The 
dependence of the single-particle crossover temperature on the interaction strength, parametrized by the exponent 
9, is given in Figure |l^. 

A point worth stressing here is that the expression (1.119) for the effective propagator is equivalent to an RPA 
approach for the interchain motion and is therefore not exact. Jsi] , |6^ , |6^ , ^ It becomes exact when interchain 
single-particle hopping has an infinite range. |]6l| Another point worth stressing is the fact that the temperature 
interval over which the single-particle dimensionality crossover takes place is not small. In effect, as already 
noticed in § 2.1, the critical ID temperature domain stretches from Ep to T making the crossover as a function 



of temperature gradual. 



5. 1 Interchain pair hopping and long-range order 

Although transverse coherent band motion is absent above T^i, a finite probability for fermions to make a jump 
on neighboring chains still survives within the thermal coherence time r ~ 1/T. When this "virtual" motion is 
combined with electron-electron interactions along the stacks, it generates an effective interchain hopping for pairs 
of particles that are created in all ID channels of correlations (Peierls, Cooper, and Landau). [l^, ^ [l5| When 
coupled with singular pair correlations along the chains, as may occur in the Cooper or Peierls channels, effective 
transverse two-particle interaction may arise and lead to the occurrence of long-range order above T^i, even when 
fermions are still confined along the chains either by thermal fluctuations or by the presence of a gap. In other 
words, interchain pair processes generate a different kind of dimensionality crossover, which in turn introduces a 
new fixed point connected to long-range order. It should be noticed, however, that long-range order is not possible 
in two dimensions if the order parameter has more than one component . ||4l[| In such a case, a finite coupling in a 
third direction is required. 



FIGURE 14. Antiferromagnetic interchain exchange V + occurring between two spins within the coherence length f . In 
the presence of a charge gap, the exchange takes place in the interval ^ — > f p ~ vp/Ap. 

To see how these new possibihties emerge from the renormalization group method described above, the quasi-lD 
propagator (1.119) is used in the evaluation of the one-loop outer shell corrections ^{Sj 2)0. c ^-^^d the perturbative 



influence of tj_ in both Peierls and Cooper channels is shown in Figure |13|. The first term is independent of t± 
and reproduces the ID one-loop corrections. The second diagram expresses the generation at each outer shell 
integration of an elementary contribution to interchain hopping of a pair of particles in the Peierls or Cooper 
channel. This outer shell generating term is given by 



Sl[r,i^]^-^z'^ii>VFY.T.f>^i^^^)0;{Q)0,iQ), (1.121) 



2 



where Q = {qo,q±),Q = {q,q±,u!m) and 



/^(Q,^) = ±2^ii^^^gl^j'cosgx. (1.122) 

The upper (lower) sign refers to the Peierls (Cooper) channel. 

This term was not present in the bare action and must be added to /is(^). Furthermore, it will be coupled to 
other terms of S giving rise to additional corrections and hence to an effective interchain two-particle action of 
the form 

s±[r.^] = E E ^'^(Q' ^)o;(Q)o^(o), (1.123) 

^' Q 

where V^(Q, i) is the pair tunneling amplitude. As a function of £, the parameter space of the action then becomes 

Hs{i) = (zGp, zt±, ziQi, 2:252, 2:353, V^,.. .) 

Among the most interesting pair tunneling processes in the Peierls channel, +,r, corresponds to an interchain 
kinetic exchange and promote site SDW or AF ordering (Figure ^4|), while V^-^^ is an interchain BOW coupling 

which is involved in spin-Peierls ordering. Both favor the transverse propagation of order at Qo = {2kp,Tr), 
which actually coincides with the best nesting vector of the whole warped Fermi surface. In the Cooper channel, 
V^P describes the interchain Josephson coupling for SS and TS Cooper pair hopping. According to Figure ^l|, 
these terms can only couple to singular intrachain correlations when Umklapp scattering is irrelevant. Long-range 
order in the Cooper channel is uniform and occurs at the wave vector Qo = (0, 0). 



FIGURE 15. Recursion relation for the interchain pair hopping term in the transverse RPA treatment (Eqn. (1.124)). 



Long-range order will occur at the critical point determined by lie — ^'0-{Ef /Tcf^) where V^(Qoi-^c) — > oo. 
A reasonable way to obtain an approximate Tc.^ is a transverse RPA approach (Fig |l^), which consists of 
treating interchain coupling in the mean-field approximation while the intrachain correlations are taken into 
account rigorously using (l.llS ) or the exact asymptotic form of the pair susceptibility when it is known. Thus 
substituting Sj ^ Sj + S± in ( 1.74 ), one g ets th e one-loop outer-shell contribution \{S'j 2 + 2S'/^2'S'±,2 + S\ 2)0- 
The first term leads to the generating term (1.121 ) to second order in zti_, whereas the second term introduces pair 



vertex corrections to the composite field O^*'' which can then be expressed in terms of given in (1.113). Finally, 
the third term is the RPA contribution that enables the free propagation of pairs in the transverse direction and 



hence the possibility of long-range order at finite temperature (Fig. 15). In all cases other than the generating 
term, the outer-shell integration in the RPA approach is made by neglecting t±_ , with the result 



^^^(Qi ^) — //J I ►A'^^ii-^AM o^''^ 



(1.124) 



for the flow equation of the interchain pair hopping. The solution of this inhomogeneous equation at I has the 
simple pole structure 





where 



(1.125) 



[^^liunyXiQ.nx.iny 



(1.126) 



Critical temperature in the random phase approximation 



Owing to the structure oiVfj,, a general analytic solution for T^^^ is impossible. However, an approximate expression 
for the critical temperature can be readily obtained, by first neglecting transients in the power laws of ^ e'"'*'^, 
z ~ e~^^, and the couplings gp.{t), and by setting to unity the term in square brackets in the integrand of (1.126). 
Thus, far from Tc^^ where the RPA term can be neglected, one can write 



F,,(Qo,^)«2 



(5^ 
El 



.((2-2^?-7^)-i(e^(2-2«-'^-)-l), 



(1.127) 



where (5^) is an average of the coupling constants over the interval [0,^]. In the presence of a gap in the charge 
or spin sector, the above scaling form at £ lv=p,a = \vl(Ef jTa^p) is governed by the strong coupling condition 

26* — 7p > 0. In strong coupling when 



26** 



7 



< , while for 



=p,(T 

^ (.u=p,<n the couplings are weak and 2 



^ I 2 — 29* — 7* |:» 1, the expression ( 1.127 ) actually coincides with so that its substitution in ( 1.125 ) leads to 
a non-self-consistent pole at 



where 



(gg)^^(C) tl 



(1.128) 



(1.129) 
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is an effective interchain pair hopping term at the step This form is reminiscent of a perturbative expansion 
in terms of the small parameter zi^/A^ indicative of interchain hopping of bound pairs. p5| , |6^ , |67t The positive 
constant 

S^. = l + g^^{^./EoV-x^.{^u) < 1 (1.130) 



corresponds to the integral of (1.126) over the interval [0,£w] at high energy. These results in the presence of a 
gap allow us to write 



J. 



(1.131) 



which is compatible with (1.27) from the extended scaling ansatz. T^2 is thus a temperature scale for the di- 
mensionality crossover for the correlation of pairs of particles induced by g±fj,, and 0j.2 = 7* is the two-particle 
exponent. |15| Tj;2 may also be equated with the onset of critical (Gaussian) fluctuations. A point that is worth 
noticing here concerns the variation of the critical temperature with the gap A^. Actually, for most interesting 
situations, the flow of coupling constants will cross the Luther-Emery line and the use of the expression (111) 
with 7* = 1 will lead to the dependence Tc^^ « z'^{£,^)t\/ A^, in which Tc.p. increases as the amplitude of the gap 
decreases (Figure |l^). In the absence of a gap, the interchain pair hopping term is still present and it can lead to 
long-range ordering. Sufficiently far from the critical point and for £{2 — 26* — 7^) ^ 1, one flnds 



^^(Qo,^) 



-2.g±^(2- 



,20-2 



(1.132) 



where in weak coupling g±f^ 
critical point will occur at 



{gl)t\/El. Substituting this expression for in the RPA expression (|1.125D , the 



EF{g^,/[l,{2-2e-^,)])-^^-^K 



(1.133) 



Therefore, in weak coupling, T^^^ 
0j,2 — 2(1 — 9) which is twice (p^i 

^2\ " 



--^ Tj,2 ~ (ff_L/i) "=^1 and one can define the two-particle crossover exponent 
(plus anomalous corrections coming from the g^). The critical temperature 
Tc,ti ~ (5^)^ z't± will then decrease when interaction decreases. This result is meaningful as long as Tr^2 > Tj.i, 
that is when the two-particle dimensionality crossover that marks the onset of Gaussian critical fluctuations still 
falls in a temperature region where the transverse single particle motion is incoherent. The interval of coupling 
constants where this occurs is relatively narrow, however, since when there is no gap t± is a relevant parameter 
and the calculated T^i increases with decreasing interactions. When the strong and weak coupling regimes are 
looked at together, a maximum of Tc,^ is therefore expected at the boundary (Figure p^ . 

Finally, it should be stressed that whenever the strong coupling condition 2 — 26* — 7 < is achieved in the 
absence of any gaps as in the case for the Tomanaga-Luttinger model, there is no maximum of Ten. Instead, it 
is found to increase monotonically with the interaction strength (see the dotted line of Figure |l6|) . B^, Q 



Response functions 

The calculation of the response functions in the presence of interchain pair hopping can be done following the 
procedure given in § 4.4. At the RPA level for the pair hopping, one must add the outer-shell contribution 
{ShS±,2)Q c to the scaling equation for Xfi in the ID regime (see Eqn. (1.113)), to get 



|lnx,(Q,^) = |lnx.(^)-^.. 



Using (1.125), the integration of this equation leads to 



Xm(^) 



(1.134) 



(1.135) 



which presents a singularity that is quadratic in (T — Tc,^)^^ at Q = Qo. The exponent 7 (cf. Eqn. ( 1.25| )) for 
the auxiliary susceptibility is therefore twice the value expected for Xp, in the transverse RPA approach. For the 
cases of two-particle crossovers in § 5.1, this expression at Qo can also be written in an extended scaling form 
similar to (1.26) 

MQo.e) « X^.iQo,e)X{Bg^^/T^^-) (1.136) 



FIGURE 16. Schematic diagram of the temperature (T) and interaction strength (9) dependence of the single-particle and 
pair dimensionality crossovers from the Luttinger liquid (LL) or Luther- Emery liquid (T < A,,) towards a Fermi liquid 
(FL) or long-range ordering (LRO). The dashed (dotted) line corresponds to the onset of LRO from a Fermi liquid (gap less 
LL in strong coupling). 



where X{y) — 1/(1 — y)^ is the RPA crossover scahng function for the auxihary susceptibihty. 

In the present transverse RPA scheme, the small longitudinal q dependence can be incorporated in through 
the substitution (1.116), while the q± dependence will come from an expansion of V"^ in ( 1.126| ) around the 
transverse ordering wavector qj_o in the channel considered. The Q dependence of will take the form 



(1.137) 



where ^0 vp/T^ is the longitudinal coherence length; the transverse coherence length ^_lo < 1 is smaller than 
the interchain distance indicating strong anisotropy. The growth of the correlation length is governed by the usual 

Gaussian expression ^ =^0^"^/^, which according to (1.3C), leads to the correlation exponent i> = 1/2. Similar 



agreement with scaling is found for the response function, for which the above RPA treatment leads to the familiar 
results j — 1 and rj = 0. 



5.2 Long-range order in the deconfined region 
Critical temperatures in the ladder approximation 

When the calculated value of T^2 becomes smaller than the single-particle deconfinement temperature T^i , the two- 
particle crossover is irrelevant and the expression (1.133) for Tc^^ becomes incorrect. In effect, thermal fluctuations 
have energies smaller than the characteristic energy related to the warping of the Fermi surface, which then 
becomes coherent. The t± expansion of Fig. ^ stops being convergent and the flow equation (1.124) must be 
modifled accordingly. Now the effective low-energy spectrum 



i?p(k) = ep(fc) — 2zt± cosfcj 



(1.138) 



is characterized by the electron- hole symmetry (nesting) £^+(k + Qo) = — £^_(k) and logarithmic corrections 
will also be present below Eo{£^i) (Figure |l^a). Similarly in the Cooper channel, the inversion property of the 
spectrum i?+(k) — i?_(— k) also leads to a logarithmic singularity in the same energy range (Figure |l|-b). WcwiU 



follow a simple two-cutoff scaling scheme, [68, 6^, |l5| in which the flows of (Eqn. (|l.l24 )) and (Eqn. ( |l.95 )) 
are stopped at i^i and their values are used as boundary conditions for the outer shell integration below E(){£^i) 
which is done with respect to constant warped energy surfaces (Figure [l7| ). The outer shell integration will be 
done at the one-loop level in the so-called ladder approximation meaning that the interference between the Cooper 
and Peierls channels will be neglected. This decoupling between the two channels at T <C T^i is actually not a 
bad approximation for Tc^^ when the couplings are attractive and favor superconductivity or when the nesting 
properties of the warped Fermi surface are good and a density-wave instability is expected. The approximation 
becomes less justified, however, if for example nesting deviations are present in the spectrum and the instability 
of the Peierls channel is weakened. In this case, it has been shown that the residual interference between the two 



FIGURE 17. (a) particle-hole pairing in the Peierls channel at Qo = {2kF,n). (b) Particle-particle (or hole-hole) pairing 
in the Cooper channel. Constant energy surfaces at ±^i5o(^ > £x^) from the Fermi surface are shown as dashed lines. 



channels gives rise to non-trivial influence on pairing which becomes non-uniform along the Fermi surface — e.g. 



unconventional superconductivity for repulsive interactions and anisotropy of the SDW gap. 1 69 
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In the ladder approximation, the evaluation of ^(S'^ 2)o,c a-t Qo in each channel allows us to write 



(1.139) 



This expression is incomplete, however, because ■^{Sj 2)0.0 evaluated at Qo, gives a logarithmic outer shell cor- 



rection of the form ^g'^di, which leads to the flow equation 



di 



^ - -(a )^ 



(1.140) 



In the action, the corresponding combination of couplings is proportional to 0*0^^ and will then add to 



Therefore it is useful to define Q^ = V^ — for which one obtains the flow equation of Fig. |l^ 



di 



The relevant parameter space for the action for Eq{£) < Eo{£^i) in the ladder approximation becomes 

/i5(^>40 = (G°(k,c^„),g^). 



The solution of (1.141) has the following simple pole structure 



1 + 5^^(40(^-40' 



(1.141) 



(1.142) 



(1.143) 



where Gfi{(x^) is the boundary condition deduced from the ID scahng equations ( |1.124|) and ( |l.95| ) at £^1. The 
presence of an instability of the Fermi liquid towards the condensation of pairs in the channel /i occurs at the 
critical temperature 

T,,^ = exp(-2/ I g^(40 I)- (1-144) 

This BCS type of instability occurs when the effective coupling constant is attractive {G^ < 0) . The phase diagram 
coincides with the one found in ID when there is perfect nesting (see Fig. pl| ). The critical temperature in each 
sector is associated with those correlations having the largest amplitude. The variation of the critical temperature 
as a function of the interaction strength is shown by the dashed line in Figure |l^. We also notice that the BCS 
expression (1.144) satisfies the extended scaling hypothesis (1.2£), namely Tc/^ oc T^i cx tY'^"^ ■ 



Response functions in the deconfined region 

By adding the coupling of fermion pairs to source fields, the calculation of response functions in the channel 
/i below T^i can readily be done. In the ladder approximation, the outer shell contraction ^{ShSj ^2)0,01 which 



FIGURE 18. Recursion relation in the ladder approximation for the effective coupling C/p = — (crossed square) in 
the channel fi below the temperature of the single-particle dimensionality crossover T^i . 



amounts making the substitution Q^^ in Fig. |l^, leads at once to the flow equation for the pair vertex part 



Using the ladder result (1.143), one gets for the auxiliary susceptibility ~ (-^^u)^ at ^ > l^i 



Xm(Qo,^) 



[1 + ^^^(40(^-40]' 



(1.146) 



which presents a power law singularity that is quadratic in [T — Tc^^)~^ as the te mpera ture approaches Tc,^ from 
above. The integration of this expression following the definition given in Eqn. ( 1.109 ) causes no difhculties and 
one finds for the response function 



Xf^iQaJ) = Xm(Qo:40 - (ttwf) 



-1 Xm(Qo^ 40(^-40 
1 + 1^^(40(^-40' 



(1.147) 



The singularity appears in the second term and is linear in {T—Tc^^) ^ , which is typical of a ladder approximation. 



The critical exponent defined in Eqn. ( 1.25 ) is 7 = 1 and corresponds to the prediction of mea n-fiel d theory. As 
for the amplitude of the singular behavior, it satisfies the extended scaling constraint shown in (1.2S). It is worth 
noting that these scaling properties are also satisfied by the auxiliary susceptibility. 

In going beyond the ladder approximation, the RG technique for fermions becomes less reliable in the treatment 
of fluctuation effects beyond the Gaussian level which are usually found sufficiently close the critical point — at 
least in three dimensions. Standard approaches of critical phenomena, such as the statistical mechanics of the 
order parameter from a Ginzburg-Landau free energy functional, can be used sufficiently close to Tc,^. 



6 Kohn-Luttinger mechanism in quasi-one-dimensional metals 

6. 1 Generation of interchain pairing channels 



The outer shell partial integration described in the RPA treatment of § 5.1 for T > T^i, namely in the presence 
of interchain pair hopping, was restricted to the channel /i corresponding to the most singular correlations of the 
ID problem (Figure |ll|). However, looking more closely at the possible contractions of fermion fields in S±^2, a 
new pairing possibility emerges as shown in Figure |l^. This consists of putting the two outer-shell fermion fields 
on separate chains i and j, which can then be equated with interchain pairing. The attraction between the two 
particles in the pair is mediated by the combination of pair hopping and intrachain correlations. Thus interchain 
pairing is similar to the Kohn-Luttinger mechanism for superconductivity by which Friedel charge oscillations 
can lead to an effective attraction between two electrons and in turn to an instability of the 3D metal at a finite 
but extremely low temperature. ||3^ In quasi-lD systems, however, the K-L mechanism has a dual nature given 
that interchain pair hopping can exist not only in the Peierls channel but also in the Cooper channel. Therefore 
for a Josephson coupling, the interchain outer-shell decomposition of Figure |l9| leads to an attraction between a 
particle and hole on different chains. In this case the attraction is mediated by the exchange of Cooper pairs. 

Coupling constants 

In the following we will restrict our description of the onset of interchain pairing correlations within the KW 
RG to the incommensurate case where only gi and 92 scattering processes are present. Thus in the temperature 



FIGURE 19. Possible outer-shell contractions for interchain pair hopping: intrachain channel ^ (left) and interchain 
channel /i (right). The black square corresponds to the pair hopping amplitude V^ij between chains i and j. 



domain above T^i, which corresponds to ^ < l^i, we have seen in § 5.1 that the outershell contraction in which 
the fermion lines are diagonal in the chain indices yields the RPA equation (1.124). As pointed out above, another 
possible contraction would be to put the two outer shell fermion fields on distinct chains. 

This is best illustrated by first rewriting the pair hopping term (1.123) in terms of chains indices. Then putting 
two fields in the outershell, we consider the sum of two contributions 



Si 



S'j 



■KVp Z 



E E ^M,.. W(o;,.(9)OM,j(g) + c.c 



-■nvp z 



E E u,..,A(■){ol^Jm^.Aq:) 



+ c.c 



(1.148) 



the first of which leads to (1.124) at the one- loop level, while the second is connected to the interchain channel 
denoted p,. The interchain composite fields are defined by 



(1.149) 



for the singlet (ISS: — 0) a-nd triplet (ITS: /ic = 1,2,3) interchain Cooper channel /xc, and by 

Opp,.,(g) = {T/LY/^ Y.^-..Uk - q)<7f^^+jj,,ik) (1.150) 



for the charge-density-wave (ICDW: fip = 0) and spin-density-wave (ISDW: = 1,2,3) pair fields of the 
interchain Peierls channel ftp. This introduces the combinations of couplings 



(1.151) 



The constants are obtained from 



MC(P)=0,ij 



1 1 



ISS (ICDW) 
ITS (ISDW). 



2 tJ-p(c)-!^^,i-j • 2 

Therefore at the one-loop level in addition to the RPA term, one obtains the term 



(1.152) 



(1.153) 



which is also logarithmic. Owing to the dependence on the chain indices of the last composite field on the right 
hand side of this equation, this term cannot be rewritten in terms of the intrachain channels and therefore stands 
as a new coupling whic h we w ill write dS±. It must be added to the parameter space of the action. To do so, it 
is convenient to recast ( 1.153 ) in terms of interchain backward and forward scattering amplitudes: 



6S± = -—TTVpZ^^ie) E E {f2,tj{£)SaiaJa 



fl.ij (^)f^ai as 



{a},q i,j 
X 



+ 9)C,a. j(fc2 - g)V'-,a3j(fc2)^+,a4,»(fcl)- (1-154) 



FIGURE 20. Perpendicular backward (gi") and and forward ((72") scattering terms generated by interchain pair tunneling. 



The outer-shell backscattering amplitude /i^^ then corresponds to intrachain momentum transfer near 2kp 
between the particles and a change in the chain indices during the process, whereas the forward generating 
amplitude f2^ij corresponds to small intrachain momentum with no change in the chain index (Figure H). The 
explicit expressions are 

2/2,»,-/i,., = ^[3F^'c/o,»,+^^'c=o.„] + ^[3^Mc#o,^,+^MP=o,.,]- (1-155) 

Since the pair hopping amplitudes are themselves generated as a function of it is clear that the f s will be 
very small except when the approach the strong coupling regime. Successive KW transformations will lead to 
the renormalization of the above interchain backward and forward scattering amplitudes which we denote Qxij 
and 92 12 ' respectively. These define the part of the action denoted S\_ , which when expressed in terms of uniform 
charge and spin fields introduced in ( |l.84 ), takes the rotationally invariant form 



s^[r,i']i = -7TVFZ-\e) J2 E ( - W) PpAq)p-pA~Q) - 3m. W ■ ^-pA-<i) )(i-156) 



p-.q y 

This expression clearly shows that interchain pair hopping generates a coupling between uniform charge and spin 
excitations of different chains. The fiow equations of these couplings will now be obtained from 5 (5^ 2)0, c ^-t the 
one-loop level. Aside from the scale dependent generating terms, the contributions consist of the same type of 
interfering Cooper and Peierls diagrams found in the purely one-dimensional problem (Fig. |^). These outer shell 
logarithmic contributions can be evaluated at once resulting in the following fiow equations for the interchain 
spin and charge variables 

|25^..-56, - yhj-ft:r (1-157) 

It is interesting to note that, following the example of the purely one-dimensional couplings, the interference 
between the interchain and Cooper channels causes these equations to be decoupled. Hence the spin and charge 
degrees of freedom on different chains are not mixed through and §2 ■ 

The solution for the charge coupling is straightforward and one obtains a parquet-type solution for the spin 
part 



where 



ffM,(^)= / dl'h^,,{f)[\^gi,Af)]\ (1.159) 



According to ( [l.l55[ ), ,fi.ij{£') and therefore gi^ij{£) are negative in the repulsive sector {gi > and 2g2 — gi > 0) 
where SDW correlations are singular (Figure |ll|) ; g^^^j {£) has then the possibility to fiow to the strong attractive 
coupling sector indicating the existence of a spin gap. Given the place held by gj'^^j- (£) in the hierarchy of generated 
couplings of the RG transformation, however, this can only occur if the the SDW pair tunneling or exchange term 
V^p^o,ij first reaches the strong coupling domain — a possibility realized as the system enters the Gaussian critical 



domain near Tx2. Since the above fiow equations (1.157) hold for all pairs of chains ij for £ < the amplitude of 
the spin gap goes to zero as the separation |i — j | ^ 00. It is worth noting that the spin gap for nearest-neighbor 
chains is reminiscent of the one occurring in the two-chain problem. Jzol [7l| 
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Interchain response functions 

To establish the nature of correlations introduced by interchain pairing, it is convenient to express outershell 
decomposition of Sj_ + S±^ in the corresponding channel 



1 



<3,M='= *J 

where we have introduced the combinations of couplings 



± 



(1.160) 



5lTS±,4j 



.9l,y - 92,ij ± UlTS,ij 



for symmetric (+) and antisymmetric (— ) interchain superconductivity and their analog expressions 



3lCDW±,ii 
5lSDW±,ij 



92,ij 



± C^ISDW,y 

for interchain density-wave. The corresponding pair fields are given by 

1 

2' 



(1.161) 



(1.162) 



(1.163) 



for symmetric (+) and antisymmetric (— ) composite fields. 

We add the coupling to source fields Sh = J2 ^*p,± ij^p,^,ij + ^^'^ one-loop corrections to the pair vertex 
part are obtained from {{S±,2 + S±^2)Sh)o.c (which is equivalent to making the substitution 5^ —> 5j^g± in 
Fig. 10), and yields the flow equations pair vertex renormalization factor Zjj^± ^j and hence that of the auxiliary 



susceptibility x^i 



^I±,^J (with z 



1): 



(1.164) 



Obtaining an expression for Xfi^.ij hi closed form for arbitrary chains i and j is not possible. However, in the case 
where intrachain channels have developed only short-range order and C/^,y is only sizable for nearest-neighbour 
chains, one can extract the dominant trend of interchain auxiliary suceptibilities in the various sectors of t he ffi.0 2 
plane. Thus by taking in Fourier space y^(Q,^) w (— l)''fMc | V^(^) | cosqx for the pair tunneling term ( 1.125 ), 
one gets V^^^j{t) = ^(-1)^^'"= | V^^it) \ where | ¥^,{1) \ = \ V;,(Qo,^) | is the maximum amplitude of V^(Q,^). This 
approximation holds sufficiently far from T^^^, namely for \V^x^^lx^l not too close to unity. The nearest-neighbour 
ij auxiliary susceptibility then becomes 



Xf.±a]{(-) 




< I V,{t) I M' 



(1.165) 



where the values of the channel index ^ and the exponent are given in (1.151-1.152). The transient amplitude 
is given by 



Xf,^ij{i) = exp ■ 



(1.166) 



where for superconductivity /3iss 



-9l.i] 



92 Ai and /3its 



9i.i] 



g2iji while for density-wave /3icdw = 



g^ij — 'igiij and /3isdw = 92 ij - The onset of interchain correlations, even though they are small, will be essentially 
determined by the most singular correlations of the intrachain channels /i (Figure pT| ) and the exponent cj^. At 
finite temperature, we can write 



X/i±,^J(7') 



(r)n [^-lix.iT))- 



(1.167) 



FIGURE 21. Phase diagram of the quasi-one-dimensional electron gas with interchain hopping and no Umklapp scattering. 
Dominant correlations of intrachain (interchain) pairing channels are given. 



In the repulsive sector where gi > 0, 2g2 — gi > 0, SDW correlations are singular (region I, Figure El|) and lead 
to an enhancement of symmetric (+) interchain singlet superconducting (188+) correlations with c^p^o ~ "^Z-^- 
should be noted here that, in this sector, symmetric interchain triplet (ITS"*") correlations are also enhanced by 
SDW corr elations with c|;^^° = 1/2. The pairing on different chains results from the exchange of spin fluctuations. 
On the T-L line at gi = and 172 > 0, SDW and CDW correlations are equally singular leading to an enhancement 
of both ISS+ and ITS+ correlations with a net exponent of unity. 



In region II of Figure 21, where gi > 0, 2^2 — 5i < 0, the TS channel is the most singularly enhanced leading 
to interchain charge-density- wave (ICDW~) correlations with cJ^^^q = 3/2. Here the pairing results from a three- 
component (triplet) Cooper pair exchange between the chains. In the same region, an enhancement, although 
smaller, of antisymmetric interchain spin-density- wave correlations (ISDW~) is also present with c^^^q ~ ^Z^- 
the T-L line at 32 < 0, however, SS and TS are equally singular so both contribute to give the same enhancement 
for ISDW~ and ICDW~. If one moves to the attractive sector (region III, Figure |2^), where gi < 0, 2^2 — 51 < 
and SS is the most singular channel along the chains, one finds that symmetric ICDW''" and antisymmetric ISDW^ 
channels have similar enhancements c|^p=o ~ '^'^10=0 ~ ^1"^- Finally in the region IV, where g\ < 0, 2g2 — 5i > 
and CDW correlations are the most singular, symmetric ITS^ correlations are enhanced in the interchain triplet 
channel with the exponent cjjp^o ~ ^f'^- 

6.2 Possibility of long-range order in the interchain pairing channels 

If the initial parameters of the model are such that the system can reach the energy or temperature domain T^i 
without developing long-range order, the pair-hopping amplitude V^i is still in the weak coupling regime. We have 



seen in §5.2, however, that infrared singularities in both the Peierls and Cooper channels persist above ^3,1 so 
that long-range order in the intrachain channel /i can occur if t/^^ < 0. This will be so in the Peierls channel when 
deviations from perfect nesting of the warped Fermi surface do not exceed a given threshold (Fig. |l^a), while 
in the Cooper channel if a magnetic field is below a critical value. Therefore sufficiently large nesting deviations 
or magnetic field must be introduced in the model in order for long-range order in the interchain channels to 
be possible. We will now illustrate how such a possibility can be realized in the repulsive case when nesting 
deviations are introduced in the model. We will thus focus our attention in the region I of the phase diagram 
where symmetric interchain singlet superconductivity (ISS+) can be stabilized (Fig. ^). A similar description 
holds for the Peierls interchain channel in the presence of a sufficiently strong magnetic field (regions II and III 
in Fig. |l|). 

Frustration of nesting is commonly introduced by considering the possibility of particles hopping to second 
nearest-neighbor chains, i.e. by adding the term — 2ij_2 cos 2A;j_ to the spectrum at ^ = 0, where t±2 is the 
amplitude of the hopping {t±2 ^ t±). The relevant parameter space of the action at the crossover value £^,1 
becomes fJ-si^x^) = {Gp{K^^n),G^L,g^±), where 

G°(k,w„) = — (1.168) 

' iujn- ep{k) + 2t*j_cosk± + 2t*j_.^cos2k± ^ ' 
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is the effective propagator in which z is the quasi-particlc weight at and ^^(2) — ^^±(2)- If one cranks up 
the ampUtude of t*j_2, nesting of the Fermi surface at Qo — {2kp,iT) is quickly suppressed so that the Tc^sdw in 



the region I of the phase diagram decreases rapidly and eventually vanishes above the critical tj 



As for the 



infrared singularity of the Cooper channel, however, it is still present and can then lead to an instability of the 
normal state. 

We will adopt here a simple ladder RG description of this instability, which neglects the interference between 
the — weakened — Peierls channel and the Cooper one. This approximation should be qualitatively valid for 
^*L2 ^ ^*f2 more quantitative description including the interference effect can be found in Ref. ^). We will 
further assume for simplicity that the enhancement of above t^'^ is negligible and the transverse short-range 
order of the chains is entirely contained in the ij dependence of 5j^g+ ^ in Eqns. (1.161) and ( 1.152| ), which is 
restricted to nearest- neighbor chains. In this way, the repulsive coupling is entirely local and does not compete 
with the interchain part: fermions on different chains are said to avoid intrachain repulsion. 134, Rn| In transverse 



Fourier space, the outer shell decomposition (1.160) in the channel ISS"*" becomes 



S 



s 



1 



ISS^ 



-TTVpZ 



,fc' 



(1.169) 



q fcj 



in which we have defined Oss(<7) — Oss(<7, k±) and where rik± = cosk± is the Fourier factor that leads to a 
global antisymmetric order parameter for interchain singlet superconductivity. The combination of couplings at £ 



(1.170) 



is obtained from the amplitude of the transverse couplings at £. In the ladder approximation, the outer shell 
integration ^((5x_2 +'5^. 2)^)0, c is done in this channel. At zero external pair momentum, this leads to the flow 
equation 

^ . 1 . _ 

(1.171) 



^ .-L l/± ^2 



which is easily integrated to give the simple pole expression 



on 

9lSS+ (^a;! ) 



(1.172) 



4^133+ 



The critical temperature for the superconducting instability is then Tc — T^i exp[— 4/gig_j_^^i )] where gf^g+ (4i ) > 
corresponds to a net attraction in the interchain superconducting channel. From ( |l.l69| ), the singlet order pa- 
rameter is of the form 

A(fci) = Ao77fc^, (1.173) 

which has nodes on the Fermi surface. 

The calculation of the interchain superconducting response causes no particular difhculty. In effect, by adding 
the coupling to a source field {Sh = k±^rjk^ hp_+ + c.c), one obtains the flow equation for the auxiliary 

susceptibility 

;lnXi33+ = ^5iis+- (1.174) 



di 



Using (1.172), this is integrated to give the double pole singularity 

Xl33+(4i) 



X133+ (^) 



(1.175) 



where the boundary condition Xiss+i^x^) is given by the solution of ( |1.164 ) at i^i. The loop integration of the 
above expression following the definition ( 1.109| ), leads to the response function 



X133+W = Xi33+(4i) - (27rwF) ^- 



Xi33+(4i)(^-4i) 

~ 4-9l33+ (4i)(^ ~ 4i) 



(1.176) 



which has the expected simple pole structure of the ladder approximation. 
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7 Summary and concluding remarks 

We have reexamined the scahng tools that enter in the description of quasi-one-dimensional interacting fermion 
systems. On phenomenological grounds, we have emphasized that the scaling ansatz can provide a qualitative 
understanding of the existence of power law behavior of various quantities in one dimension while its extension 
for non zero interchain coupling has allowed us to describe the different forms of dimensionality crossover in 
the quasi-one-dimensional case. A key feature of this phenomenology is the existence of fluctuations of single 
fermion and pair degrees of freedom on many length or energy scales, the origin of which can be microscopically 
understood using the renormalization group language. We found it useful to reconsider the formulation of this 
methodology and we have seen that the adoption of a classical Wilson procedure for the Kadanoff transformation 
has a molding influence on the RG flow beyond the one-loop level. Thus the constraint of a sharp natural cutoff ko 
on the fermion spectrum which keeps the loop momenta in the outer energy shells inevitably leads to a non-local 
flow of the parameters that define the action. Non-locality results from the inclusion of many-particle marginal 
interactions in the outer shell integration. These are not present in the bare action but are generated along the RG 
flow. The resulting structure of outer shell integration finds an interesting parallel with the functional description 
of the renormalization group equations proposed by Polchinski in the framework of the theory, and used 
by Zanchi and Schulz in the context of two-dimensional interacting fermion systems. ||20| 

It is only when the marginally relevant or irrelevant couplings are treated in the local approximation that 
purely logarithmic scaling is recovered. Although this has only been verified at the two-loop loop level and in 
the one-dimensional case, similar conclusions are expected to hold at higher order and in higher dimensions as 
well. The procedure described in this review casts light on a previous alternative formulation of the RG which 
is known to offer a more direct route to logarithmic scaling at high order but which on the other hand requires 
the relaxation of some of the standard rules of the Wilson method. It also highlights how the difficulties linked 
to non-locality of the flow equations are essentially hidden in older versions of the renormalization group. 

The RG method is widely recognized as a tool well suited to the study of the influence of interchain hopping on 
the stability of the Luttinger and Luther-Emery liquid states in one dimension. We have reviewed the conditions 
under which the coupling between stacks can make these non-Fermi liquid states unstable. A possible instabil- 
ity occurs when the Luttinger liquid yields to the formation of a Fermi liquid. The recovery of quasi-particle 
excitations, here called a dimensionality crossover for the fermion quantum coherence, squares well with the pic- 
ture provided by the extended scaling ansatz. A second instability emerges when the kinematics of interchain 
single-particle hopping is combined with intrachain pair correlations. Effective pair hopping processes between 
chains then form and their coupling to singular correlations along the chains lead to a distinct (two-particle) 
dimensionality crossover. This is physically meaningful as long as the corresponding temperature scale is above 
the one predicted for single particles. Such a situation is favored by cranking up interactions which on the one 
hand increases pair correlations and on the other hand reduces single particle coherence. The temperature scale 
for two-particle dimensionality crossover marks the onset of long-range order which accords well with the ex- 
tended scaling ansatz. For sufficiently weak interactions, however, single-particle deconfinement occurs first and 
the possibility for long-range order remains as an instability of the Fermi liquid. In this case, the ladder level of 
the RG gives a BCS character to the transition. 

The RG method can also be useful in studying how the Kohn-Luttinger mechanism for unconventional pairing 
works in the quasi-one-dimensional case. We have seen that this mechanism opens up channels of correlation 
which pair particles on different stacks as a result of pair hopping processes between chains. Since these are 
present in both the Pcicrls and Cooper channels, the Kohn-Luttinger mechanism has a dual nature. Thus it was 
shown not only that the exchange of density-wave correlations can lead to interchain Cooper pairing but also 
that interchain Peierls pairing between a particle and a hole is also possible when both particles exchange Cooper 
pairs. Long-range order in these unconventional channels can only be achieved from the Fermi liquid state in 
which enough nesting deviations or magnetic field are present to frustrate the formation of long-range order in 
the the primary channels of correlation along the chains. 

The results presented in this review by no means exhaust the range of possibilities supplied by the RG method. 
We will end this article by giving two specific examples. The first concerns the issue of single-particle deconfinement 
that leads to the recovery of a Fermi liquid. This question is of great interest especially in connection with the 
normal state properties obeserved in concrete realizations of quasi-one-dimensional materials like the organic 
conductors, [p^ Q We have seen that the critical domain of the primary (Luttinger liquid) fixed point in a quasi- 
one-dimensional system essentially covers the whole temperature range of interest. This favors in turn a slow and 
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gradual crossover to a Fermi liquid, which means in practice that transitory effects are likely to be important in 
the description of these materials. Although crossover scaling effects have been underestimated in the treatment 
given above, they are not beyond the reach of the RG method. 

Another promising avenue is the application of the RG method to interacting fermions on a lattice. The 
Kadanoff transformation can actually be generalized to fermions with the full tight binding spectrum. Transients 
to the continuum limit can then be considered in a quantitative way. In the case of the fermion representation of 
the one-dimensional spin chain, for example, this has been shown to be an essential ingredient in the description 
of the thermodynamics. 
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1 One-particle self-energy at the two-loop level 

1.1 Backward and forward scattering contributions 

Following the partial trace operation in the absence of umklapp scattering, three-particle interactions are generated 
from the outer shell contraction SS'"^^ = \{S'j i)a.c- The corresponding expression reads 

X V'-.a!,(fc2)V'-,a3(fc2)'0+,a4(fcl) 

X V+,a^(fc'l)V'+,a4(fcl)V'-,a3(fc2) (177) 
where A^^_^,j(^j_3) = g{a}(^j-3).9{a'}(^i-3)'5aia^, ^{a,a'}i^J-3) = 5{a} (^j-3)5{a'} (^j-s)'^^^^^ , = HEo/Eo{ej- 

and 5{q} is given by (1.62). Here the summation ^ ~ covers all uJm but q is restricted such that ki -\- q (^2 — q) 

i n th e propagator G^j. (G*^) is in the outer momentum shell at the step (j — 2)d£ for j > 3. The two terms in 
( |l77|) , which depends on the step j, correspond to the diagrams of Figure ^a. At the next partial integration 
{j — l)d£ , the above contributions will in turn be contracted leading to a second fermion line at the adjacent 
momentum shell (Figure H-b). When this set of contractions is repeated as a function of j, contributions add to 
give an effective two-particle interaction of the form 

N-l 

6S',^ E('^^i!'2)o.c 

= (--i^)V5IE E E E^{«,.'}(V3){'5.^..GO(fc + 5)G"(fc'-g) 

k,q' {a,a'} f., g 

X (fc + ? + q)f.,o.', O^'-J^ 9)V^-,"3 (fc')V'+,«. (fc) 

+ 6^,^, G°(fc-l-§)G°(fc'-g') 

(fc)} 
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+ i^-^r^ E EEE E^Ko'}(V3)'5.;.„,G'[(fci+^G^(fci) 



N-l 



{a,a.'} 



3=3 



fci 



(178) 



for A'*' and an analogous expression for A^ . In the first two terms, the summation £p covers all fermion frequencies 
LOn' and its momentum interval is such that fc' — q is in the outer shell while k' is still in the inner shell. These 
expressions correspond to the two series of diagrams shown in Figure ||-c, while the last term (not shown in 
Figure ^c) only involves outgoing fermions near —kp. The final contraction of the cascade is obtained by putting 
k' in the outer momentum shell at NdH. 



A,2/0,c 



rp2 N-l 

^i^^Ffj^Y. E E E (9!{ej-3) + g^iij-s) - 51(^,-3)52(^,-3)) 



{ [g° (fc + q)G°_{k' - q)G°_{k') + Gl{k + q)G"_{k' + q)G"_{k')] Jfc)V+,a(fc) 
Glik' + q)G°_(k - q)G°_(k') + ...] C«(fc)V'-,a(fc)}. 



(179) 



The first two terms in the sum lead to a self-energy correction for fermions (on the +kp branch) the momenta, 
k, of which are located at the top of the inner shell. The outer shell contractions leading to the last terms are 
self-energy corrections for fermions on the —kp branch (only one such term is shown explicitly). Owing to the 
outer shell constraints put on the momentum, however, when the summations over fermion and boson frequencies 
are carried out in the low temperature limit, all these —kp branch terms give a vanishingly small contribution in 
comparison with the first two terms and can be safely ignored. 



('5'5'^,2)o,, 



J^Eo(In) ^ J-^En(lN)-^Eo(t^-3) 

X {gl{l,-z) + 9l{l,-z) - 31(^,-3)52(^,-3)) 



^Eo(tN-l) ^ J ^Eo{lN) + ^Eo(lj-2) 

^ r+o.{k)Mk)- 



[G\{k)]-^-2vpq 




The integrand turns out to be independent of e' = t-{k') and we will define SEo/2 — J_{ 



-^EoiiN) 



n 



^Eo{£n-i) 
Eo{In) 



de' 



de' . Now assuming that the coupling constants vary slowly as a function of j, this expression can 



be evaluated in the local approximation (Figure 




-Eoit)) 



+ 



i{Eo{i)+Eo) 



^iEo{l)+Eo) 
Eo(l) 



-^(5?(^) + 52'(^)-5l(%2(^))E 



{gliQ + gli^,) - gii^M^,)) 

+ / 



d{vpq) 



[GO(fc)]- 
d[vpq) 



[G\{k)]- 



2vpq 



(gl (^) + 52' (^) - 51 (^)52 (i) W)de E [G° (fc)] - 1 V4a (fc) V'+a (fc) 



^;jfc)V+a(fc) 

zvpq 

r+o.m+aik) 

(181) 



to leading order in [G^(fc)]-i. Here £g ^ IhEq/ \vpq\, t{£) = 1 - 2/(1 + e/) and £'o(£iv-3) ^ Eo{£) in the limit 
d£ —> 0. The final expression then becomes the outer shell contribution of the renormalization coefficient z~^{£) 
of [G!|_]~^ for the inner shell state k in the local approximation. The presence of the factor t{£) indicates that the 
above result is not logarithmic over the entire range, especially at small £ where transients or scaling deviations 
exist, [p^ These become vanishingly small at large £ and are neglected to logarithmic accuracy. A similar cascade 
of contractions applies to the second term of (177) with A~ and this leads to the same renormalization for the 
one-particle propagator near —kp. 



1.2 Umklapp contribution 

When umklapp scattering is present, the contraction SSu^ 
three-particle interaction at the step {j — 2)d£ with j > 3: 



/0,c: 



that involves g^ alone, generates a 
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{k,Q'} Q' 

X V-,ai(^l)V'+,Q'i(^l)'0+,a^(fc2) 

+ i^^Ff^ EE A^(^,_3) G^(^i+g) Ca;(fc-Q')C,a,(fcl+Q + 0'-G)^;,a.(fc2-Q + G) 
{fc,Q} Q 

X V+,a;(fci)^-,ai(A:i)V'-,«.(fe) (182) 

where A^(^j_3) = .g|(^j_3), Q = (2/ci? + g, i^m) and G = (4fcF,0). Here £ ;^ ( £ ?;,) consists in summing all 
boson frequencies and those q [q') values of the momentum transfer for which fci + 9 (fc^ — q') is in the outershell 
at step (j — 2)M. Following the example of the incommensurate case, the next two contractions of the first term 



of (182) will lead to corrections to the one-particle term at +fcF with internal lines in three successive outer shells 



and the external k located at the top of the inner shell: 



N 2 ^ 

Y.i5S'!^)o.. = (vrt;^)2^EE E E.93(^.-3)G"+(fc' + Q)G^(P)G'L(fc + Q-G)^;„(fc)V'+„(fc) (183) 



L 

J=3 a 3=3 J., ^ 



The explicit evaluation runs parallel to the incommensurate case in (17E-180). Thus in the local approximation, 
the variation of gs{£j^3) with j is neglected and the value of the coupling is taken at {N — l)d£ — > £. The outer 
shell umklapp contribution to the single particle term at the inner shell state k of the action then becomes 

^ 1 

E (<^^>o,c ^ -9l{lW)MY.[Gl{k)]-'rUky^+c,ik). (184) 

i=3 a 



This result, which adds to z ^{d£), corresponds to the diagram of Figure ^e. Using the second term of (182) a 
similar expression holds for fermions at —kp. 
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